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Results of searches for heavy stable charged particles produced in pp collisions at 
y/s = 7 and 8 TeV are presented corresponding to an integrated luminosity of 5.0 fb _1 
and 18.8 fb -1 , respectively. Data collected with the CMS detector are used to study the 
momentum, energy deposition, and time-of-flight of signal candidates. Leptons with 
an electric charge between e/3 and 8e, as well as bound states that can undergo charge 
exchange with the detector material, are studied. Analysis results are presented for 
various combinations of signatures in the inner tracker only, inner tracker and muon 
detector, and muon detector only. Detector signatures utilized are long time-of-flight 
to the outer muon system and anomalously high (or low) energy deposition in the 
inner tracker. The data are consistent with the expected background, and upper limits 
are set on the production cross section of long-lived gluinos, scalar top quarks, and 
scalar x leptons, as well as pair produced long-lived leptons. Corresponding lower 
mass limits, ranging up to 1322 GeV/c 2 for gluinos, are the most stringent to date. 
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1 Introduction 

Many extensions of the standard model (SM) include heavy, long-lived, charged particles that 
have speed v significantly less than the speed of light c HH3 or charge Q not equal to the 
elementary positive or negative charge ±le [4-8], or both. With lifetimes greater than a few 
nanoseconds, these particles can travel distances comparable to the size of modern detectors 
and thus appear to be stable. These particles, generically referred to as heavy stable charged 
particles (HSCP), can be singly charged (|Q| = le), fractionally charged (|Q| < le), or multiply 
charged (|Q| > le). Without dedicated searches, HSCPs may be misidentified or even com- 
pletely missed, as particle identification algorithms at hadron collider experiments generally 
assume signatures appropriate for SM particles, e.g., v ~ c and Q = or ±le. Additionally, 
some HSCPs may combine with SM particles to form composite objects. Interactions of these 
composite objects with the detector may change their constituents and possibly their electric 
charge, further limiting the ability of standard algorithms to identify them. 

For HSCP masses greater than 100GeV/c 2 , a significant fraction of particles produced at the 
Large Hadron Collider (LHC) have /3 (= v / c) values less than 0.9. These HSCPs can be iden- 
tified by their longer time-of-flight (TOF) to outer detectors or their anomalous energy loss 
(dE/dx). The dE/dx of a particle depends on both its electric charge (varying as Q 2 ) and its j6. 
The dependence of dE/dx on these variables is described by the Bethe-Bloch formula j9|. This 
dependence can be seen in Fig. [IJ which shows a dE/dx estimate versus momentum for tracks 
from data and Monte Carlo (MC) simulations of HSCP signals with various charges. In the 
momentum range of interest (10-1000 GeV/c), SM charged particles have a relatively flat ion- 
ization energy loss and /3 values very close to one. to massive particles with |Q| = le, particles 
with |Q| > le, and low-momentum particles with |Q| < le. On the other hand, searching for 
candidates with lower dE/dx yields sensitivity to high-momentum particles with |Q| < le. 

Previous collider searches for HSCPs have been performed at LEP lH0Hl3l , HERA HH, the 
Tevatron IT5H18I , and the LHC IH9TI261 . The results from such searches have placed important 
bounds on beyond the standard model (BSM) theories [27, 28 j, such as lower limits on the 
mass of gluinos, scalar top quarks (stops), and pair-produced scalar t leptons (staus) at 1098, 
737, and 223GeV/c 2 , respectively. Presented here are several searches for singly, fractionally, 
and multiply charged HSCPs using data collected with the Compact Muon Solenoid (CMS) 
detector during the 2011 {y/s = 7TeV, 5.0 fbr 1 ) and 2012 (y/s = 8TeV, 18.8 fb _1 ) data taking 
periods. 

2 Signal benchmarks 

The searches presented here are sensitive to a wide variety of signals of new charged mas- 
sive particles. Several BSM models are used to benchmark the sensitivity. The HSCPs can be 
classified as either lepton-like or hadron-like. Lepton-like HSCPs interact primarily through the 
electromagnetic force, while hadron-like HSCPs additionally interact through the strong force 
and form bound states with SM quarks (or gluons) called R-hadrons [29 j. The R-hadrons can 
be charged or neutral. Strong interactions between the SM quarks and detector material in- 
crease energy loss and can lead to charge exchange, e.g., conversion of charged R-hadrons into 
neutral ones (and vice-versa). There is some uncertainty in the modeling of R-hadrons' strong 
interactions with detector material. For this analysis, two separate models are considered: (1) 
the model described in Refs. 11301 l3Tf , referred to as the cloud model, and (2) a model in which 
any strong interaction results in a neutral R-hadron [32], referred to as the charge-suppressed 
model. The cloud model envisions the R-hadron as composed of a spectator HSCP surrounded 
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2 Signal benchmarks 




Figure 1: Distribution of 7/„ a dE/dx estimator that is defined in Section 3.1 versus particle 
momentum for y/s = 8 TeV data (left) and also including MC simulated HSCP candidates of 
different charges (right). Tracks with 2.8 < Ij, < 3.0 MeV/cm are excluded by preselection 
requirements, as discussed in Section |4j 

by a cloud of colored, light constituents. The charge-suppressed model results in essentially 
all R-hadrons being neutral by the time they enter the muon system. For each of the models 
considered, particle interactions with the CMS apparatus and detector response are simulated 
using Geant4 v9.2 11331 |34H . To produce the effect of multiple interactions per bunch crossing 
(pileup), simulated minimum bias events are overlaid with the primary collision. 

The minimal gauge-mediated supersymmetry breaking (GMSB) model [35| predicts the grav- 
itino to be the lightest supersymmetric particle (LSP) and allows for the next-to-lightest su- 
persymmetric particle (NLSP) to be long-lived because of the weakness of the gravitational 
coupling, which governs the decay of the NLSP to the LSP. For this analysis the NLSP is taken 
to be a lepton-like stau (ti) with an assumed lifetime that exceeds the time-of-flight through the 
CMS detector. For y/s = 7 (8) TeV simulation, pythia v6.422 (v6.426) ESI is used to model both 
Drell-Yan production of a Ti pair (direct pair-production) and production of heavier supersym- 
metric particles whose decay chains lead to indirect stau production. Events with T\ masses 
in the range 100-494 GeV/c 2 are generated using line 7 of the "Snowmass Points and Slopes" 
benchmarks |37|. They correspond to N = 3 chiral SU(5) multiplets added to the theory at 
a scale F from 60 to 320 TeV [35 1 depending on the % mass, and an effective supersymmetry- 
breaking scale of A = F / 2. All points have a value of 10 for the ratio of neutral Higgs field 
vacuum expectation values (tan /3), a positive sign for the Higgs-Higgsino mass parameter 
(sgn(^)), and a value of 10 4 for the ratio of the gravitino mass to the value it would have if 
the only supersymmetry-breaking scale were that in the messenger sector (Cg rav )- The parti- 
cle mass spectrum and the decay table were produced with the program ISASUGRA version 
7.69 [38]. Theoretical production cross sections ((%) for staus are calculated at next-to-leading 
order (NLO) with PROSPINO v2.1 H39| . 

R-hadron signals from gluino (g) and scalar top (ti) pair production are studied using PYTHIA 
v6.442 (v8.153) H36l |40| for y/s = 7 (8) TeV generation. Stop pair production is modeled for 
masses in the range 100-1000 GeV/c 2 . For g production, split supersymmetry fl4Tl |42| is mod- 
eled by setting the squark masses to greater than 10 TeV/ c 2 and generating g masses of 300- 
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1500GeV/c 2 . The fraction / of gluinos hadronizing into g-gluon bound states is unknown. 
These neutral states would not leave a track in the inner detectors. Therefore, several scenar- 
ios are considered for the singly charged analysis: / = 0.1, 0.5, and 1.0. In the extreme case 
where / = 1.0, R-hadrons are always neutral in the inner tracker, but a fraction of them may 
interact with the detector material and be electrically charged during their passage through the 
muon system. Gluino and scalar top pair production cross sections are calculated at NLO plus 
next-to-leading logarithmic (NLL) accuracy with PROSPINO v2.1 Il43l450l . 

The last of the signal samples studied is the modified Drell-Yan production of long-lived lep- 
tons. In this scenario, new massive spin-1/2 particles may have an electric charge different 
than \Q\ = le and are neutral under SU(3)c and SU(2)l) therefore they couple only to the 
photon and the Z boson via U(l) couplings [51J. Signal samples are simulated using PYTHIA 
v6.422 (v6.426) |36| for y/s = 7 (8) TeV. The analysis uses simulations of |Q| = e/3 and 2e/3 for 
masses of 100-600 GeV/c 2 , of \Q\ ranging from le to 6e for masses of 100-1000 GeV/c 2 , and of 
| Q| = 7e and 8e for masses of 200-1000 GeV/c 2 . 

The CTEQ6L1 parton distribution functions (PDF) [52] are used for the sample generation. 



The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
interaction point, the x axis pointing to the center of the LHC ring, the y axis pointing up 
(perpendicular to the plane of the LHC ring), and the z axis along the counterclockwise-beam 
direction. The polar angle 9 is measured from the positive z axis and the azimuthal angle <p in 
the x-y plane. The pseudorapidity is given by rj = — ln[tan(0/2)]. 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame- 
ter. Within the field volume are a silicon pixel and strip tracker, a lead tungstate crystal electro- 
magnetic calorimeter, and a brass and scintillator hadron calorimeter. Muons are measured in 
gas-ionization detectors embedded in the steel flux-return yoke of the magnet. Extensive for- 
ward calorimetry complements the coverage provided by the barrel and endcap detectors. The 
inner tracker measures charged particles within the pseudorapidity range \rj\ < 2.5. It consists 
of 1440 silicon pixel and 15 148 silicon strip detector modules and is located in the 3.8 T field of 
the solenoid. The inner tracker provides a transverse momentum (pj) resolution of about 1.5% 
for 100 GeV/c particles. Muons are measured in the pseudorapidity range \tj \ < 2.4, with detec- 
tion planes made using three technologies: drift tubes (DT), cathode strip chambers (CSC), and 
resistive plate chambers (RPC). The muon system extends out to eleven meters from the inter- 
action point in the z direction and seven meters radially. Matching tracks in the muon system 
to tracks measured in the silicon tracker results in a transverse momentum resolution between 
1 and 5%, for pj values up to 1 TeV/c. The first level (LI) of the CMS trigger system, composed 
of custom hardware processors, uses information from the calorimeters and muon detectors 
to select events of interest. The high level trigger (HLT) processor farm further decreases the 
event rate from around 100 kHz to around 300 Hz for data storage. A more detailed description 
of the CMS detector can be found in Ref. [53 J. 

3.1 The dE/dx measurements 

As in Ref. |19| , dE/dx for a track is estimated as: 
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3 The CMS detector 



where N is the number of measurements in the silicon-strip detectors and c, is the energy loss 
per unit path length in the sensitive part of the silicon detector of the z'th measurement; If, has 
units MeV/cm. In addition, two modified versions of the Smirnov-Cramer-von Mises |54ll55| 
discriminator, I as (I' as ), are used to separate SM particles from candidates with large (small) 
dE/dx. The discriminator is given by: 



ias - N x I i2N + L 



2N 



(2) 



where P, (P-) is the probability for a minimum ionizing particle (MIP) to produce a charge 
smaller (larger) or equal to that of the z'th measurement for the observed path length in the 

detector, and the sum is over the measurements ordered in terms of increasing P>'' . 

As in Ref . [19J, the mass of a \Q\ = le candidate particle is calculated based on the relationship: 

m 

I h = K^ + C, (3) 

where the empirical parameters K = 2.559 ± 0.001 MeV • c 2 /cm and C = 2.772 ± 0.001 MeV/cm 
are determined from data using a sample of low-momentum protons in a minimum-bias dataset. 

The number of silicon-strip measurements associated with a track, 15 on average, is sufficient 
to ensure good dE/dx and mass resolutions. 

3.2 Time-of-flight measurements 

As in Ref. [24[, the time-of-flight to the muon system can be used to discriminate between 
/3 « 1 particles and slower candidates. The measured time difference (St) of a hit relative to 
that expected for a /3 = 1 particle can be used to determine the particle 1 / /3 via the equation: 

cS 

1/0 = 1 + -j-, (4) 

where L is the flight distance from the interaction point. The track 1 / j5 value is calculated as 
the weighted average of the 1 / /3 measurements from the DT and CSC hits associated with the 
track. A description of how the DT and CSC systems measure the time of hits is given below. 

As tubes in consecutive layers of DT chambers are staggered by half a tube, a typical track 
passes alternatively to the left and to the right of the sensitive wires in consecutive layers. The 
position of hits is inferred from the drift time of the ionization electrons assuming the hits come 
from a prompt muon. For a late arriving HSCP, the delay will result in a longer drift time being 
attributed, so hits drifting left will be to the right of their true position while hits drifting right 
will be to the left. The DT measurement of St then comes from the residuals of a straight line fit 
to the track hits in the chamber. The weight for the z'th DT measurement is given by: 

Wi _ (^i^L, (5) 



DT 



where n is the number of (p projection measurements found in the chamber from which the 
measurement comes and o~dt = 3 ns is the time resolution of the DT measurements. The factor 
(n — 2)/n accounts for the fact that residuals are computed using two parameters of a straight 
line determined from the same n measurements (the minimum number of hits in a DT chamber 
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needed for a residual calculation is n = 3). Particles passing through the DTs have on average 
16 time measurements. 

The CSC measurement of St is found by measuring the arrival time of the signals from both 
the cathode strips and anode wires with respect to the time expected for prompt muons. The 
weight for the z'th CSC measurement is given by: 

L? 

m = (6) 

where a - ,, the measured time resolution, is 7.0 ns for cathode strip measurements and 8.6 ns 
for anode wire measurements. Particles passing through the CSCs have on average 30 time 
measurements, where cathode strip and anode wire measurements are counted separately. 

The uncertainty {(T\/p) on 1/0 for the track is found via the equation: 



f (1/fr- 1/0)2 XWj 



N - 1 



where 1/0 is the average 1/0 of the track and N is the number of measurements associated 
with the track. 

Several factors including the intrinsic time resolution of the subsystems, the typical number 
of measurements per track, and the distance from the interaction point lead to a resolution of 
about 0.065 for 1/0 in both the DT and CSC subsystems over the full rj range. 



4 Data selection 

Multiple search strategies are used to separate signal from background depending on the na- 
ture of the HSCP under investigation. 

• For singly charged HSCPs, 

o the "tracker+TOF" analysis requires tracks to be reconstructed in the in- 
ner tracker and the muon system, 

o the "tracker-only" analysis only requires tracks to be reconstructed in the 
inner tracker, and 

o the "muon-only" analysis only requires tracks to be reconstructed in the 
muon system. 

• For fractionally charged HSCPs, the "fractionally charged" (|Q| < le) analysis only 
requires tracks to be reconstructed in the inner tracker and to have a dE/dx smaller 
than SM particles. 

• For multiply charged HSCPs, the "multiply charged" (|Q| > le) analysis requires 
tracks to be reconstructed in the inner tracker and the muon system. The analysis is 
optimized for much larger ionization in the detector compared to the tracker+TOF 
analysis. 

The tracker-only and muon-only cases allow for the possibility of charge flipping (charged to 
neutral or vice versa) within the calorimeter or tracker material. The muon-only analysis is the 
first CMS search that does not require an HSCP to be charged in the inner tracker. The singly, 
multiply, and fractionally charged analyses feature different selections, background estimates, 
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4 Data selection 



and systematic uncertainties. The preselection requirements for the analyses are described 
below. 

All events must pass a trigger requiring either the reconstruction of (i) a muon with high pj or 
(ii) large missing transverse energy (E™ 1SS ) defined as the magnitude of the vectorial sum of the 
transverse momenta of all particles reconstructed by an online particle-flow algorithm l[56| at 
theHLT. 

The LI muon trigger allows for late arriving particles (such as slow moving HSCPs) by accept- 
ing tracks that produce signals in the RPCs within either the 25 ns time window corresponding 
to the interaction bunch crossing or the following 25 ns time window. For the data used in 
this analysis, the second 25 ns time window is empty of proton-proton collisions because of the 
50 ns LHC bunch spacing during the 2011 and 2012 operation. 

Triggering on E™ ss allows for some recovery of events with hadron-like HSCPs in which none 
of the R-hadrons in the event are charged in both the inner tracker and the muon system. The 
E™ 1SS in the event arises because the particle-flow algorithm rejects tracks not consistent with a 
SM particle. This rejection includes tracks reconstructed only in the inner tracker with a track 
pi much greater than the matched energy deposited in the calorimeter [57] as would be the case 
for R-hadrons becoming neutral in the calorimeter, and tracks reconstructed only in the muon 
system as would be the case for R-hadrons that are initially neutral. Thus, only the energy these 
HSCPs deposit in the calorimeter, roughly 10-20 GeV, will be included in the E™ 1SS calculation. 
In events in which no HSCPs are reconstructed as muon candidates, significant E™ 1SS will result 
if the vector sum of the HSCPs' momenta is large. The E™ ss trigger will collect these events, 
allowing for sensitivity to HSCP without a muon-like signature. 

For all the analyses, the muon trigger requires pj > 40GeV/c measured in the inner tracker 
and the E™ 1SS trigger requires E2p iss > 150 GeV at the HIT. The muon-only analysis uses the 
same two triggers, and a second muon trigger that requires both pj > 70 GeV/ c (using only 
the muon system) and E™ 1SS > 55 GeV. For the first part of the 2012 data (corresponding to 
an integrated luminosity of 700 pb -1 ), the requirement was E2p iss > 65 GeV. Using multiple 
triggers in all of the analyses allows for increased sensitivity to HSCP candidates that arrive 
late in the muon system and to hadron-like HSCPs that are sometimes charged in only one of 
the inner tracker and muon system and sometimes charged in both. The muon-only analysis 
uses only y/s = 8 TeV data as the necessary triggers were not available in 2011. 

For the tracker-only analysis, all events are required to have a track candidate in the region 
\t] | < 2.1 with pj > 45GeV/c (as measured in the inner tracker). In addition, a relative un- 
certainty in pi (o~p T / pi) less than 0.25 and a track fit x 2 P er number of degrees of freedom 
(rid) less than 5 is required. Furthermore, the magnitudes of the longitudinal (d z ) and trans- 
verse (d X y) impact parameters are both required to be less than 0.5 cm. The impact parameters 
d z and d xy are both defined with respect to the primary vertex that yields the smallest \d z \ 
for the candidate track. The requirements on the impact parameters are very loose compared 
with the resolutions for tracks in the inner tracker. Candidates must pass isolation require- 
ments in the tracker and calorimeter. The tracker isolation requirement is Epx < 50GeV/c 
where the sum is over all tracks (except the candidate) within a cone about the candidate track 
AR = \/ (Ar/) 2 + (Acp) 2 < 0.3 radians. The calorimeter isolation requirement is E/p < 0.3 
where E is the sum of energy deposited in the calorimeter towers within AR < 0.3 and p is 
the candidate track momentum reconstructed from the inner tracker. Candidates must have 
at least two measurements in the silicon pixel detector and at least eight measurements in the 
combination of the silicon strip and pixel detectors. In addition, there must be measurements 
in at least 80% of the silicon layers between the first and last measurements on the track. To 
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reduce the rate of contamination from clusters with large energy deposition due to overlap- 
ping tracks, a "cleaning" procedure is applied to remove clusters in the silicon strip tracker 
that are not consistent with passage of only one charged particle (e.g., a narrow cluster with 
most of the energy deposited in one or two strips). After cluster cleaning, there must be at least 
six measurements in the silicon strip detector that are used for the dE/dx calculation. Finally, 
J/j > 3MeV/cm is required. 

The tracker+TOF analysis applies the same criteria, but additionally requires a reconstructed 
muon matched to a track in the inner detectors. At least eight independent time measurements 
are needed for the TOF computation. Finally, 1//3 > 1 and Ci/a < 0.07 are required. 

The muon-only analysis uses separate criteria that include requiring a reconstructed track in 
the muon system with pj > 80GeV/c within \t]\ < 2.1. The relative resolution in pj is ap- 
proximately 10% in the barrel region and approaches 30% for \t]\ > 1.8 ||58| . However, charge 
flipping by R-hadrons can lead to an overestimate of pj. The effect is more pronounced for 
gluinos, where all of the electric charge comes from SM quarks. The measured curvature in 
the muon system for gluinos is 60-70% smaller than would be expected for a muon with the 
same transverse momentum. The candidate track must have measurements in two or more 
DT or CSC stations, and \d z \ and \d xy \ < 15cm (calculated using tracks from the muon system 
and measured relative to the nominal beam spot position rather than to the reconstructed ver- 
tex). The requirements on \d z \ and \d xy \ are approximately 90% and 95% efficient for prompt 
tracks, respectively. To reject cosmic ray muons, candidates are removed if there is a muon 
segment both with rj within ±0.1 of — ?/ can d, where rj can d is the pseudorapidity of the HSCP 
candidate, and with \5(f>\ > 0.3 radians, where Sep is the difference in cp between the candi- 
date and the muon segment. Also candidates compatible with vertically downward cosmic 
ray muons, 1.2 < \cp\ < 1.9 radians, are rejected. To reject muons from adjacent beam crossings, 
tracks are removed if their time leaving the interaction point as measured by the muon system 
is within ±5 ns of a different LHC beam crossing. This veto makes the background from muons 
from such crossings negligible while removing very little signal. Finally, the same quality re- 
quirements used in the tracker+TOF analysis are applied in the muon-only 1 / /3 measurement. 

The fractionally charged search uses the same preselection criteria as the tracker-only analysis 
except that J/, is required to be <2.8MeV/cm. An additional veto on cosmic ray muons rejects 
candidates if a track with pj > 40 GeV/ c is found on the opposite side of the detector (AR > 

7T-0.3). 

The multiply charged particle search uses the same preselection as the tracker+TOF analysis 
except that the E/p selection is removed. Furthermore, given that a multiply charged parti- 
cle might have a cluster shape different from that of a |Q| = le particle, the cluster cleaning 
procedure is not applied for the multiply charged analysis. 

The preselection criteria applied on the inner tracker track for the analyses are summarized in 
Table [l] while the criteria on the muon system track are summarized in Table |2j 

5 Background prediction 

Candidates passing the preselection criteria (Section [4]) are subject to two (or three) additional 
selection criteria to further improve the signal-to-background ratio. For all of the analyses, 
results are based upon a comparison of the number of candidates passing the final section 
criteria with the number of background events estimated from the numbers of events that fail 
combinations of the criteria. 
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5 Background prediction 



Table 1: Preselection criteria on the inner tracker track used in the various analyses as defined 
in the text. 





IQI <le 


tracker+TOF tracker-only Q > \e 


\n\ 


<2.1 


p T (GeV/c) 


>45 


d z and d xy (cm) 


<0.5 


^Pt/Pt 


<0.25 


Track x l l n & 


<5 


# Pixel hits 


>1 


# Tracker hits 


>7 


Frac. Valid hits 


>0.8 


Zpf(AR < 0.3) (GeV/c) 


<50 


# dE/dx measurements 


>5 


dE/dx strip shape test 


yes 


no 


E cal (AR <0.3)/p 


<0.3 




I h (MeV/cm) 


<2.8 


>3.0 


AR to another track 


<7t-0.3 





Table 2: Preselection criteria on the muon system track used in the various analyses as defined 
in the text. 





tracker+TOF \Q\ > le 


muon-only 


# TOF measurements 


>7 




<0.07 


1/(6 


>1 


n 




<2.1 


p T (GeV/c) 




>80 


d z and d xy (cm) 




<15 


# DT or CSC stations 




>1 


Opp. segment rj\ difference 




>0.1 






<1.2 0R >1.9 


\St\ to another beam crossing (ns) 




>5 
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Tracker + TOF CMS is = 8 TeV, L = 1 8.8 fb' 



Tracker + TOF CMS tfs = 8 TeV, L= 18.8 fb 1 
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Figure 2: Measured I as (left) and 1/B (middle) distributions for several p T ranges and measured 
1/B distributions for several I as ranges (right). Results are for the tr acker +TOF selection at 
^Js = 8 TeV. The lack of variation of the distributions for different ranges of the other variables 
demonstrates the lack of strong correlation between 1/B, I as , and pj. 

The background expectation in the signal region, D, is estimated as D = BC/ A, where B (C) 
is the number of candidates that fail the first (second) criteria but pass the other one and A 
is the number of candidates that fail both criteria. The method works if the probability for 
a background candidate to pass one of the criteria is not strongly correlated with whether it 
passes the other criteria. The lack of strong correlation between the selection criteria is evident 
m Fig.|2j All tracks passing the preselection enter either the signal region D or one of the control 
regions that is used for the background prediction. 

For the tracker-only analysis, the two chosen criteria are pj and I as . Threshold values (pj > 
70GeV/c and I as > 0.4) are fixed such that failing candidates passing only pj (I as ) fall into 
the B (C) regions. The B (C) candidates are then used to form a binned probability density 
function in 4(p) such that, using the mass determination (Eq. ([3])), the full mass spectrum of 
the background in the signal region D can be predicted. The rj distribution of candidates at low 
dE/dx differs from the distribution of the candidates at high dE/dx. To correct for this effect, 
events in the C region are weighted such that the rj distribution matches that in the B region. 

For the tracker+TOF analysis, three criteria are used, pj, l as , and 1/ B, creating eight regions la- 
beled A through H. The final threshold values are selected to be px > 70 GeV/c, I as > 0.125, and 
1/ B > 1.225. Region D represents the signal region, with events passing all three criteria. The 
candidates in the A, F, and G regions pass only the 1/B, l as , and pj criteria, respectively, while 
the candidates in the B, C, and H regions fail only the pj, I as , and 1/B criteria, respectively. The 
E region fails all three criteria. The background estimate can be made from several different 
combinations of these regions. The combination D = AGE /E 2 is used because it yields the 
smallest statistical uncertainty. Similar to the tracker-only analysis, events in the G region are 
reweighted to match the t] distribution in the B region. From a consideration of the observed 
spread in background estimates from the other combinations, a 20% systematic uncertainty is 
assigned to the background estimate. The 20% systematic uncertainty is also assigned to the 
background estimate for the tracker-only analysis. 

In order to check the background prediction, loose selection samples, which would be domi- 
nated by background tracks, are used for the tracker-only and tracker+TOF analyses. The loose 
selection sample for the tracker-only analysis is defined as pj > 50GeV/c and I as > 0.10. The 
loose selection sample for the tracker+TOF analysis is defined as pj > 50 GeV/c, I as > 0.05, and 
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5 Background prediction 
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Figure 3: Observed and predicted mass spectra for candidates entering the tracker-only (left 
column) or tracker+TOF (right column) signal region for the loose selection. The expected dis- 
tribution for a representative signal scaled to the integrated luminosity is shown as the shaded 
histogram. The top (bottom) row is for ifs = 7 (8) TeV. 

l/j8 > 1.05. Figure [3] shows the observed and predicted mass spectra for these samples. 

The muon-only analysis uses the pi and 1/j8 criteria for the ABCD method. The final selections 
are px > 230GeV/c and 1/jS > 1.4. It has been found that these variables are correlated with 
\rj | and with the number of muon stations used to fit the candidate. Therefore, the background 
prediction is performed in six separate bins (2/3/4 muon stations in central (\t]\ < 0.9) and 
forward (0.9 < \t]\ < 2.1) regions). The final result is computed from a sum of these six 
bins. The systematic uncertainty in this background estimate is determined by defining four 
additional regions A',B',C, and D' . Events in B' (A') pass (fail) the pj requirement, but with 
0.8 < 1//3 < 1.0, while events in D' (C') pass (fail) the pj requirement with 1//3 < 0.8. Two 
complementary predictions now become possible, D = CB' / A' and D = CD'/C'. From a 
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Figure 4: Observed and predicted numbers of tracks in both the control region with 1/ B < 1 
(left) and the signal region (right) as functions of the 1/ B threshold and for two different px 
thresholds for the muon-only analysis at \/s = 8 TeV. Only statistical uncertainties are shown. 

consideration of the spread of the three estimates, a systematic uncertainty of 20% is assigned 
to the background estimate for the muon-only analysis using this method. 

The muon-only analysis also has background contributions from cosmic ray muons even af- 
ter the previously mentioned cosmic ray muon veto requirements are applied. The number of 
cosmic ray muons expected to pass the selection criteria is determined by using the sideband 
region of 70 < \d z \ < 120cm. To increase the number of cosmic ray muons in the sideband 
region, the veto requirements are not applied here. To reduce the contamination in the side- 
band region due to muons from collisions, the tracks are required to not be reconstructed in the 
inner tracker. The number of tracks (N) in the sideband with 1 / B greater than the threshold 
is counted. To determine the ratio (R fl ) of candidates in the signal region with respect to the 
sideband region, a pure cosmic ray sample is used. The sample is collected using a trigger 
requiring a track from the muon system with px > 20GeV/c, rejecting events within ±50 ns of 
a beam crossing and events triggered as beam halo. The cosmic ray muon contribution to the 
muon-only analysis signal region is determined as N x Rp, . A similar procedure is used to sub- 
tract the estimated cosmic ray muon contribution to the A, B, and C regions prior to estimating 
the collision muon background in the D region. The cosmic ray muon contribution to the signal 
region constitutes approximately 60% of the total expected background. The systematic uncer- 
tainty in the cosmic ray muon contribution is determined by comparing estimates using \d z \ 
ranges of 30-50 cm, 50-70 cm, 70-120 cm, and > 120 cm. It is found to be 80%. Figure [4] shows 
the numbers of predicted and observed candidates in both the control region with 1 / B < 1 and 
the signal region for various px and 1/ B thresholds for the \/s = 8 TeV data. The number of 
predicted events includes both the cosmic ray muon and collision muon contributions. Only 
statistical uncertainties are shown. 

The multiply charged analysis uses the I as and 1/ B criteria. Since the default track reconstruc- 
tion code assumes |Q| = le for px determination, the transverse momentum for |Q| > le 
particles is underestimated by a factor of le/ \Q\. Therefore px is not used in the final selection. 
In addition, while dE/dx scales as Q 2 , the dynamic range of the silicon readout of individual 
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5 Background prediction 




Figure 5: Observed and predicted numbers of tracks as a function of the I as threshold for two 
different 1 / /3 thresholds at = 7 TeV (left) and 8 TeV (right) for the multiply charged analysis. 
Only statistical uncertainties are shown. 

strips saturates for energy losses ~3 times that of a /5 ~ 1, |Q| = le particle. This results 
in (i) a bias in the reconstructed mass that is therefore not used for this analysis and (ii) a 
larger incompatibility of dE/dx measurements with respect to the MIP hypothesis, increasing 
the separation power of the dE/dx discriminator for multiply charged particles. The systematic 
uncertainty in the background estimate for the multiply charged analysis is determined by the 
same method that is used for the collision muon background in the muon-only analysis except 
with px changed to be I as . The two complementary estimates from the < 1.0 region lead 
to a 20% uncertainty. Figure [5] shows the numbers of predicted and observed candidates for 
various I as and thresholds. Only the statistical uncertainties are shown. 

The fractionally charged analysis uses the same method to estimate the background as the 
tracker-only analysis, replacing the l as variable with l' as and not applying a mass requirement. 
The systematic uncertainty in the prediction is taken from the tracker+TOF analysis. In addi- 
tion, the cosmic ray muon background is studied, since particles passing through the tracker 
not synchronized with the LHC clock often produce tracker hits with low energy readout. The 
cosmic ray muon background is found to be small and a 50% uncertainty is assigned to this pre- 
diction. The numbers of predicted and observed candidates for various pj and I as thresholds 
can be seen in Fig. [6] Only the statistical uncertainties are shown. 

For each analysis, fixed selections on the appropriate set of I as , I' as , Vi, and 1 / /3 are used to 
define the final signal region (and the regions for the background prediction). These values are 
chosen to give discovery potential over the signal mass regions of interest. For the tracker-only 
and tracker+TOF analyses, an additional requirement on the reconstructed mass is applied. 
The mass requirement depends upon the HSCP signal. For a given model and HSCP mass, 
the range is M reC o - 2cr to 2 TeV/c 2 where M reco is the average reconstructed mass for the given 
HSCP mass and a is the expected resolution. Both M rec o and a are determined from simulation. 

Table [3] lists the final selection criteria, the predicted numbers of background events, and the 
numbers of events observed in the signal region. Agreement between prediction and observa- 
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Figure 6: Observed and predicted numbers of tracks as a function of the l' as threshold for two 
different pj thresholds at \/s = 7 TeV (left) and 8 TeV (right) for the fractionally charged analy- 
sis. Only statistical uncertainties are shown. 

tion is seen over the full range of analyses. Figure [7] shows the observed and predicted mass 
distributions for the tracker-only and tracker+TOF analyses with the tight selection. The bump 
at lower mass values expected from the signal MC is due to the saturation of the strip electronic 
readout. 



6 Systematic uncertainties 

The sources of systematic uncertainty include those related to the integrated luminosity the 
background prediction, and the signal acceptance. The uncertainty in the integrated luminosity 
is 2.2% (4.4%) at y/s = 7 (8) TeV H59llgQl . The uncertainties in the background predictions are 
described in Section |U 

The signal acceptance is obtained from MC simulations of the various signals processed through 
the full detector simulation (Section |2|. Systematic uncertainties in the final results are domi- 
nated by uncertainties in the differences between the simulation and data evaluated in control 
samples. The relevant differences are discussed below. A summary of the systematic uncer- 
tainties is given in Table |4j 

The trigger acceptance is dominated by the muon triggers for all the models except for the 
charge-suppressed scenarios. The uncertainty in the muon trigger acceptance arises from sev- 
eral effects. A difference of up to 5% between data and MC simulation events has been ob- 
served II58I . For slow moving particles, the effect of the timing synchronization of the muon 
system is tested by shifting the arrival times in simulation to match the synchronization offset 
and width observed in data, resulting in an acceptance change of 2% (4%) for \Js = 7 (8) TeV. 
For the |Q| < le samples, an additional uncertainty arises from the possibility of losing hits 
because their ionization in the muon system is closer to the hit threshold. The uncertainty in 
the gains in the muon system is evaluated by shifting the gain by 25%, yielding an acceptance 
change of 15% (3%) for \Q\ = e/3(2e/3) samples. The uncertainty in the E™ ss trigger accep- 
tance is found by varying, at HLT level, the energy of simulated jets by the scale uncertainties. 
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Figure 7: Observed and predicted mass spectra for candidates entering the tracker-only (left 
column) or tracker+TOF (right column) signal region for the tight selection. The expected dis- 
tribution for a representative signal scaled to the integrated luminosity is shown as the shaded 
histogram. The top (bottom) row is for y/s = 7 (8) TeV. 
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Table 3: Results of the final selections for the predicted background and the observed numbers 
of events. The uncertainties include both statistical and systematic contributions. 





Selection criteria 


Number ( 
y/s — 7 lev 


)f events 

y/s — 8 lev 




Pt 
(GeV/c) 


jC) 




Mass 
(GeV/c 2 ) 


Pred. Obs. 


Pred. Obs. 


Tracker-only 


>70 


>0.4 




>0 
>100 
>200 
>300 
>400 


7.1 ±1.5 8 
6.0 ±1.3 7 
0.65 ±0.14 
0.11 ±0.02 
0.030 ± 0.006 


33 ± 7 41 
26 ± 5 29 
3.1 ±0.6 3 
0.55 ±0.11 1 
0.15 ±0.03 


Tracker+TOF 


>70 


>0.125 


> 1.225 


>0 
>100 
>200 
>300 


8.5 ±1.7 7 
1.0 ±0.2 3 
0.11 ±0.02 1 
0.020 ± 0.004 


44 ± 9 42 
5.6 ±1.1 7 
0.56 ±0.11 
0.090 ± 0.02 


Muon-only 


>230 




>1.40 






6±3 3 


IQI >le 




>0.500 


> 1.200 




0.15 ±0.04 


0.52 ±0.11 1 


\Q\ < le 


>125 


>0.275 






0.12 ±0.07 


1.0 ±0.2 



The E^ 1SS uncertainty for a/s = 7TeV samples is estimated to be less than 2% for all scenar- 
ios except for the charge-suppressed ones, where it is estimated to be <5%. For sfs = 8TeV 
samples it is less than 1% for all samples. 

The energy loss in the silicon tracker is important for all the analyses except for the muon-only 
one. Low-momentum protons are used to quantify the agreement between the observed and 
simulated distributions for I/, and I as - The dE/dx distributions of signal samples are varied in 
the simulation by the observed differences, in order to determine the systematic uncertainty. 
Because the fractionally charged analysis is also sensitive to changes to the number of hits on 
the track, track reconstruction is also performed after shifting dE/dx. The uncertainty in the 
signal acceptance varies by less than 24% for the |Q| = le samples, being less than 10% for 
all masses above 200GeV/c 2 . For the |Q| < le samples, the effect of the dE/dx shift and the 
track reconstruction combined is 25% (<10%) for |Q| = e/3(2e/3). The |Q| > le samples have 
sufficient separation of the signal from the final I as selection that the effect of the dE/dx shift is 
negligible. 

The Z boson decays to muons are used to test the MC simulation of the 1 / f> measurement. At 
\Js = 7TeV, the 1//3 measurement was observed to have a disagreement of 0.02 in the CSC 
system and 0.003 in the DT system. At \fs = 8 TeV a disagreement of 0.005 is observed for 
both systems. The uncertainty in the signal acceptance is estimated to be between and 15% 
by shifting 1//3 by these amounts. The uncertainty is generally less than 7% except for the 
high-charge /low-mass samples in the multiply charged analysis. 

The uncertainties in the efficiencies for muon reconstruction [58J and track reconstruction |6T| 
are less than 2% each. The track momentum uncertainty for the muon-only analysis is deter- 
mined by shifting 1/px of muon system tracks by 10%. For all other analyses, the momentum 
from the inner tracker track is varied as in Ref. |24|. The uncertainty is estimated to be <5% for 
all but the |Q| < le samples, low-mass |Q| > le samples, and the muon-only scenarios, where 
the uncertainty is less than 10%. 

The uncertainty in the number of pileup events is evaluated by varying by 5-6% the minimum 
bias cross section used to calculate the weights applied to signal events in order to reproduce 
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7 Results 



Table 4: Systematic uncertainties for the various HSCP searches. 



Signal acceptance 




Tracker-only 


Tracker+TOF 


\J\ > ±e 


Muon-only 


— Trigger acceptance 


<iO /o 


</ /o 


</ /o 


<b /o 


no/ 
/ /o 


— Track momentum scale 


< 10% 


<4% 


< 3% 


<10% 


<10% 


— Track reconstruction eff. 


<25% 


<2% 


<2% 


<2% 




— Ionization energy loss 


<18% 


<15% 


<12% 




— Time-of-flight 






<2% 


<15% 


<3% 


— Muon reconstruction eff. 






2% 


2% 


2% 


— Pile-up 


<2% 


<2% 


<2% 


<2% 


<4% 


— Detector material 


<1% 


<1% 


<1% 


20% 


<1% 


Total signal acceptance 


<31% 


<32% 


<31% 


<29% 


<13% 


Expected collision bckg. 


20% 


20% 


20% 


20% 


20% 


Expected cosmic ray bckg. 


50% 








80% 


Integrated luminosity 


2.2% (4.4%) for yjs = 7 


(8)TeV 



the pileup observed in data. This results in uncertainties due to pileup of less than 4%. 

The uncertainty in the amount of material in the detector simulation results in an uncertainty 
in the signal trigger and reconstruction acceptance, particularly for the |Q| > le samples. This 
is evaluated by increasing the amount of material in the hadronic calorimeter by a conservative 
5% [62j. Since it was not practical to evaluate the effect in detail for each value of Q considered, 
the largest change in signal acceptance observed (~20%) was assigned to all |Q| > le scenarios. 
The change in signal acceptance is < 1% for all |Q| < le scenarios. 

The total systematic uncertainty in the signal acceptance for the tracker-only analysis is less 
than 32% and is less than 11% for all of the gluino and scalar top cases. For the tracker +TOF 
analysis it is less than 15% for all cases except for |Q| = 2e/3, where the uncertainty ranges 
from 15% to 31%, being larger at low masses. The muon-only analysis has an uncertainty in 
the signal acceptance in the range of 7-13%. The multiply charged analysis has an uncertainty 
in the signal acceptance in the range of 21-29% for |Q| > le samples and 7-13% for |Q| = le 
samples with both being larger at low masses. The fractionally charged analysis has an uncer- 
tainty in the signal acceptance of 31% and 12% for \Q\ = e/3 and 2e/3 samples, respectively. 

The statistical uncertainty in the signal acceptance is small compared to the total systematic 
uncertainty for all the cases except for the low-mass highly charged scenarios, where the low 
acceptance leads to a statistical uncertainty that is comparable with the systematic uncertain- 
ties. For example, in the |Q| = 6e, M = 100GeV/c 2 signal, the statistical uncertainty is as high 
as 30%. In all cases the statistical uncertainty is taken into account when setting limits on signal 
cross sections. 

7 Results 

No significant excess of events is observed over the predicted backgrounds. The largest excess 
for any of the selections shown in Table [3]has a significance of 1.3 standard deviations. Cross 
section limits are placed at 95% confidence level (CL) for both \fs = 7 and 8 TeV using the 
CL S approach |63l 1641 where p-values are computed with a hybrid bayesian-frequentist tech- 
nique [65] that uses a lognormal model Il54l |55| for the nuisance parameters. The latter are 
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Table 5: Expected and observed cross section limits and the signal acceptance for gluino signals 
at \/s = 7 and 8 TeV, as well as the ratio of the cross section limit to the theoretical value for 
the combined dataset. The limit on the ratio for the muon-only analysis uses only \/s = 8 TeV 
data. The minimum reconstructed mass required (M req.) for each sample in the tracker-only 
analysis is also given. 



Mass 


Mreq. 


£7(pb)(v^ = 7TeV) 


cr(pb) U/s = 8TeV) 


a/a th (7+8TeV) 


(GeV/c 2 ) 


(GeV/c 2 ) 


Exp. Obs. Acc. 


Exp. Obs. Acc. 


Exp. 


Obs. 






Gluino (/ = 0.1) - 


— tracker-only analysis 






300 


>100 


0.0046 0.0063 0.17 


0.0055 0.0055 0.15 


4.0 x 10~ b 


4.6 x 10~ 5 


700 


>370 


0.0028 0.0029 0.21 


0.00081 0.00084 0.19 


0.0017 


0.0018 


1100 


>540 


0.0039 0.0040 0.15 


0.0010 0.0011 0.14 


0.098 


0.10 


1500 


>530 


0.0088 0.0092 0.066 


0.0021 0.0022 0.073 


5.1 


5.2 




Gluino charge-suppressed (/ 


= 0.1) — tracker-only analysis 




300 


>130 


0.035 0.034 0.021 


0.013 0.013 0.048 


0.00011 


0.00012 


700 


>340 


0.012 0.013 0.046 


0.0020 0.0021 0.077 


0.0044 


0.0043 


1100 


>410 


0.018 0.019 0.033 


0.0025 0.0026 0.061 


0.24 


0.25 


1500 


>340 


0.034 0.035 0.017 


0.0045 0.0045 0.035 


11 


11 






Gluino (/ = 0.5) - 


— muon-only analysis 






300 






0.0060 0.0065 0.058 


5.8 x W~ b 


6.3 x 10- 5 


700 






0.0026 0.0022 0.12 


0.0062 


0.0051 


1100 






0.0024 0.0020 0.13 


0.24 


0.20 


1500 






0.0030 0.0024 0.11 


7.5 


6.2 






Gluino (/ = 1.0) - 


— muon-only analysis 






300 






0.0070 0.0075 0.050 


6.8 x lO" 5 


7.3 x lO"* 


700 






0.0032 0.0027 0.10 


0.0075 


0.0063 


1100 






0.0030 0.0025 0.11 


0.30 


0.25 


1500 






0.0037 0.0031 0.087 


9.5 


7.9 



the integrated luminosity the signal acceptance, and the expected background in the signal 
region. The uncertainty in the theoretical cross section is not considered as a nuisance parame- 
ter. For the combined dataset, the limits are instead placed on the signal strength (ji = a/ c t h). 
Limits on the signal strength using only the 8 TeV dataset for the muon-only analysis are also 
presented. The observed limits are shown in Figs. [8 -[10 for all the analyses along with the theo- 
retical predictions. For the gluino and scalar top pair production, the theoretical cross sections 
are computed at NLO+NLL H45H48I using PROSPINO |66l with CTEQ6.6M PDFs |67|. The un- 
certainty bands on the theoretical cross sections include the PDF uncertainty, as well as the u s 
and scale uncertainties. Mass limits are obtained from the intersection of the observed limit 
and the central value of the theoretical cross section. For the combined result, the masses for 
which the signal strength is less than one are excluded. 

From the final results, 95% CL limits on the production cross section are shown in Tables |5j |6j 
|7j and [8] for gluino, scalar top, stau, and for Drell-Yan like production of fractionally, singly, or 
multiply charged particles, respectively. The limits are determined from the numbers of events 
passing all final criteria (including the mass criteria for the tracker-only and tracker+TOF anal- 
yses). Figure [8] shows the limits as a function of mass for the tracker-only and tracker+TOF 
analyses. The tracker-only analysis excludes gluino masses below 1322 and 1233GeV/c 2 for 
/ = 0.1 in the cloud interaction model and charge-suppressed model, respectively. Stop masses 
below 935 (818) GeV/c 2 are excluded for the cloud (charge-suppressed) models. In addition, the 
tracker+TOF analysis excludes r% masses below 435(339) GeV/c 2 for the direct+indirect (direct 
only) production. Drell-Yan signals with |Q| = 2e/3 and |Q| = le are excluded below 220 and 
574GeV/c 2 , respectively. 

The limits from the muon-only analysis for the scalar top and the gluino with various hadron- 
ization fractions / are shown in Fig. [9l The muon-only analysis excludes gluino masses below 
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Figure 8: Upper cross section limits at 95% CL on various signal models for the tracker-only 
analysis (left column) and tracker+TOF analysis (right column). The top row is for the data at 
y^s = 7 TeV, the middle row is for the data at ^fs = 8 TeV, the bottom row shows the ratio of 
the limit to the theoretical value for the combined dataset. In the legend, 'CS' stands for the 
charge-suppressed interaction model. 
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Muon - Only CMS = 8 TeV, L = 1 8.8 fb" 1 Muon - Only CMS = 8 TeV, L = 18.8 fb 1 




Figure 9: Upper cross section limits at 95% CL on various signal models for the muon-only 
analysis for the data at i/s = 8 TeV (left). Limits on the signal strength (ji = a I c th ) for the same 
data (right). 

Table 6: Expected and observed cross section limits and the signal acceptance for scalar top 
signals at \fs = 7 and 8 TeV, as well as the ratio of the cross section limit to the theoretical value 
for the combined dataset. The minimum reconstructed mass required (M req.) for each sample 
in the tracker-only analysis is also given. 



Mass 
(GeV/c 2 ) 


Mreq. 

(GeV/c 2 ) 


tr(pb)(v^ = 7TeV) 
Exp. Obs. Acc. 


(7(pb) (^ = 8 TeV) 
Exp. Obs. Acc. 


a/a th (7+8TeV) 
Exp. Obs. 


Stop — tracker-only analysis 


200 
500 
800 


>0 
>120 
>330 


0.0080 0.0088 0.14 
0.0024 0.0025 0.24 
0.0021 0.0022 0.28 


0.0051 0.0050 0.18 
0.0027 0.0034 0.23 
0.00072 0.00073 0.22 


0.00026 0.00029 
0.022 0.026 
0.21 0.22 


Stop charge-suppressed — tracker-only analysis 


200 
500 
800 


>0 
>120 
>270 


0.063 0.075 0.020 
0.0086 0.0089 0.066 
0.0071 0.0076 0.079 


0.018 0.026 0.050 
0.0068 0.0081 0.10 
0.0019 0.0023 0.10 


0.0011 0.0014 
0.062 0.070 
0.61 0.74 



Table 7: Expected and observed cross section limits and the signal acceptance for stau signals 
at -y/s = 7 and 8 TeV, as well as the ratio of the cross section limit to the theoretical value for the 
combined dataset. The minimum reconstructed mass required (M req.) for each sample in the 
tracker+TOF analysis is also given. 



Mass 
(GeV/c 2 ) 


M req. 

(GeV/c 2 ) 


£r(pb)(v^ = 7TeV) 
Exp. Obs. Acc. 


^(pb)(^ = 8TeV) 
Exp. Obs. Acc. 


cr/cr th (7+8 TeV) 
Exp. Obs. 


Direct +indirect produced stau — tracker+TOF analysis 


126 
308 
494 


>40 
>190 
>330 


0.0046 0.0035 0.29 
0.00094 0.0015 0.63 
0.00079 0.00084 0.74 


0.0042 0.0042 0.25 
0.00029 0.00028 0.56 
0.00023 0.00024 0.66 


0.0074 0.0065 
0.16 0.21 
1.9 1.9 


Direct produced stau — tracker+TOF analysis 


126 
308 
494 


>40 
>190 
>330 


0.0056 0.0046 0.26 
0.0011 0.0017 0.54 
0.00084 0.00088 0.69 


0.0044 0.0043 0.24 
0.00035 0.00035 0.46 
0.00025 0.00026 0.61 


0.18 0.16 
0.62 0.66 
4.7 5.0 
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Figure 10: Upper cross section limits at 95% CL on various signal models for the fractionally 
charged analysis (left column) and multiply charged analysis (right column). The top row is 
for the data at \fs = 7 TeV, the middle row is for the data at \/s = 8 TeV, the bottom row shows 
the ratio of the limit to the theoretical value for the combined dataset. 
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Table 8: Expected and observed cross section limits and the signal acceptance for the Drell-Yan 
like production of fractionally, singly, and multiply charged particles at y/s = 7 and 8TeV, as 
well as the ratio of the cross section limit to the theoretical value for the combined dataset. The 
minimum reconstructed mass required (M req.) for each sample in the tracker+TOF analysis is 
also given. 



Mass 
(GeV/c 2 ) 


M req. 

(GeV/c 2 ) 


cr(pb)(v^ = 7TeV) 
Exp. Obs. Acc. 


LT(pb)(v^ = 8TeV) 
Exp. Obs. Acc. 


cr/cT th (7+8TeV) 
Exp. Obs. 


Drell-Yan like production |Q| = e/3 — |Q| < le analysis 


100 
200 
400 




0.019 0.022 0.029 
0.0094 0.011 0.060 
0.0058 0.0066 0.098 


0.016 0.014 0.012 
0.0066 0.0058 0.030 
0.0041 0.0035 0.048 


0.19 0.17 
1.2 0.99 
15 13 


Drell-Yan like production \Q \=2e/ 3 — |Q| < le analysis 


100 
200 
400 




0.0011 0.0012 0.53 
0.00071 0.00076 0.81 
0.00083 0.00090 0.68 


0.00042 0.00039 0.45 
0.00027 0.00024 0.68 
0.00033 0.00034 0.56 


0.0015 0.0013 
0.014 0.012 
0.35 0.31 


Drell-Yan like production \ Q\ = le — tracker+TOF analysis 


200 
500 
800 


>120 
>340 
>530 


0.0015 0.0036 0.41 
0.00098 0.0010 0.60 
0.0010 0.0010 0.58 


0.00077 0.0013 0.36 
0.00028 0.00029 0.56 
0.00030 0.00031 0.52 


0.019 0.040 
0.41 0.44 
7.5 8.1 


Drell-Yan like production |Q = 2e — Q > le analysis 


200 
500 
800 




0.0016 0.0016 0.36 
0.00098 0.0010 0.59 
0.0011 0.0011 0.55 


0.00050 0.00073 0.33 
0.00029 0.00042 0.56 
0.00029 0.00042 0.56 


0.0028 0.0040 
0.11 0.15 
1.9 2.5 


Drell-Yan like production Q |= 3e — \Q > le analysis 


200 
500 
800 




0.0031 0.0034 0.18 
0.0012 0.0013 0.47 
0.0012 0.0013 0.47 


0.00090 0.0013 0.18 
0.00035 0.00051 0.46 
0.00033 0.00048 0.49 


0.0023 0.0032 
0.059 0.083 
0.99 1.4 


Drell-Yan like production Q = 4e — Q > le analysis 


200 
500 
800 




0.0082 0.0088 0.069 
0.0018 0.0020 0.31 
0.0017 0.0018 0.34 


0.0021 0.0030 0.078 
0.00051 0.00074 0.32 
0.00045 0.00064 0.37 


0.0031 0.0045 
0.048 0.068 
0.75 1.0 


Drell-Yan like production Q = 5e — Q > le analysis 


200 
500 
800 




0.030 0.032 0.019 
0.0035 0.0037 0.16 
0.0026 0.0027 0.22 


0.0066 0.0096 0.025 
0.00086 0.0013 0.19 
0.00066 0.00096 0.24 


0.0064 0.0092 
0.052 0.073 
0.71 1.0 


Drell-Yan like production Q = 6e — Q > le analysis 


200 
500 
800 




0.17 0.19 0.0032 
0.0079 0.0086 0.072 
0.0054 0.0056 0.11 


0.026 0.037 0.0064 
0.0016 0.0023 0.10 
0.0011 0.0016 0.15 


0.018 0.026 
0.055 0.077 
0.81 1.1 


Drell-Yan like production Q = 7e — Q > le analysis 


200 
500 
800 




1.1 1.2 0.00053 
0.022 0.023 0.026 
0.012 0.012 0.049 


0.086 0.13 0.0019 
0.0028 0.0042 0.057 
0.0017 0.0024 0.099 


0.047 0.063 
0.079 0.11 
0.96 1.3 


Drell-Yan like production Q = 8e — Q > le analysis 


200 
500 
800 




9.8 11 7.1 x 10~ 5 
0.068 0.072 0.0084 
0.028 0.030 0.020 


0.26 0.38 0.00064 
0.0051 0.0074 0.032 
0.0027 0.0039 0.062 


0.11 0.15 
0.11 0.16 
1.0 1.4 
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Figure 11: Lower mass limits at 95% CL for various models compared with previously pub- 
lished results |[T9T|26| . The model type is given on the x-axis (left). Mass limits are shown for 
Drell-Yan like production of fractionally, singly, and multiply charged particles (right). These 
particles were assumed to be neutral under SU(3)c and SU(2)l- 

1250(1276) GeV/c 2 for / = 1.0(0.5). 

Figure [lO] shows the limits applied to the Drell-Yan production model for both the fractionally 
charged and multiply charged analyses. The fractionally charged analysis excludes masses 
below 200 and 480GeV/c 2 for |Q| = e/3 and 2e/3, respectively. The multiply charged analysis 
excludes masses below 685, 752, 793, 796, 781, 757, and 715GeV/c 2 for \Q\ = 2e, 3e, 4e, 5e, 6e, 
7e, and 8e, respectively. The multiply charged analysis is not optimized for singly charged 
particles but can set a limit and is able to exclude masses below 517GeV/c 2 . As expected, this 
limit is not as stringent as the one set by the tracker+TOF analysis but does allow results to be 
interpolated to non-integer charge values (such as |Q| = 3e/2, Ae/3) using results from the same 
analysis. 



The mass limits for various signals and electric charges are shown in Fig. 11 and are compared 
with previously published results. 

The mass limits obtained for the reanalyzed = 7 TeV dataset are similar to the previously 
published CMS results except for the stau scenarios, where the new mass limits are slightly 
worse than the previously published ones. This result is a consequence of having a common 
selection for all mass points and models in contrast to what was done in Ref. I|24| , where the 
selection was optimized separately for each mass point and model. The mass limit for \Q\ < 
le samples are significantly improved with respect to Ref. fl23"|, thanks to a different analysis 
approach and to the use of the I' as likelihood discriminator that maximally exploits all the dE/dx 
information associated with a track. 
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8 Summary 

A search for heavy stable charged particles has been presented, based on several different sig- 
natures, using data recorded at collision energies of 7 and 8 TeV. Five complementary analy- 
ses have been performed: a search with only the inner tracker, a search with both the inner 
tracker and the muon system, a search with only the muon system, a search for low ionizing 
tracks, and a search for tracks with very large ionization energy loss. No significant excess is 
observed in any of the analyses. Limits on cross sections are presented for models with the 
production of gluinos, scalar tops, and staus, and for Drell-Yan like production of fractionally, 
singly, and multiply charged particles. The models for R-hadron-like HSCPs include a varying 
fraction of g— gluon production and two different interaction schemes leading to a variety of 
non-standard experimental signatures. Lower mass limits for these models are given. Gluino 
masses below 1322 and 1233GeV/c 2 are excluded for / = 0.1 in the cloud interaction model 
and the charge-suppressed model, respectively. For / = 0.5 (1.0), gluino masses below 1276 
(1250) GeV/c 2 are excluded. For stop production, masses below 935 (818) GeV/c 2 are excluded 
for the cloud (charge-suppressed) models. In addition, these analyses exclude Ti masses below 
435 (339)GeV/c 2 for the direct+indirect (direct only) production. Drell-Yan like signals with 
\Q\ = e/3, 2e/3, le, 2e, 3e, 4e, 5e, 6e, 7e, and 8e are excluded with masses below 200, 480, 574, 
685, 752, 793, 796, 781, 757, and 715 GeV/c 2 , respectively. These limits are the most stringent to 
date. 
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